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ABSTRACT 

We examine the importance of secular stellar mass loss for fueling ongoing star formation in disk galaxies 
during the late stages of their evolution. For a galaxy of a given stellar mass, we calculate the total mass loss 
rate of its entire stellar population using star formation histories derived from the observed evolution of the 
M„-star formation rate relation, along with the predictions of standard stellar evolution models for stellar mass 
loss for a variety of initial stellar mass functions. Our model shows that recycled gas from stellar mass loss can 
provide most or all of the fuel required to sustain the current level of star formation in late type galaxies. Stellar 
mass loss can therefore remove the tension between the low gas infall rates that are derived from observations 
and the relatively rapid star formation occurring in disk galaxies. For galaxies where cold gas infall rates have 
been estimated, we demonstrate explicitly that stellar mass loss can account for most of the deficit between 
their star formation and infall rates. 

Subject headings: cosmology: theory - galaxies: evolution - galaxies: formation - stars:formation - methods: 
numerical 



1. INTRODUCTION 

An understanding of how galaxies get their gas is key to 
a complete picture of galaxy formation and evolution. At 
all redshifts, star forming galaxies appear to be living fast 
and dangerously: they are observed to be converting their 
gas to stars at a rapid rate that would lead to exhaustion 
of their cold gas reservoi rs in about two billion years with- 
out a fresh supply of gas ( Kennicu tq|1998t iLerov et aLll2008t 
iBigiel et all 120081: iGenzel et alJl2010t iDaddi et al.ll2010bllal) . 
Observations, however, also show that star forming galaxies 
are common at all redshifts, and theoretical models indicate 
that star formation rates of the most actively star f orming sys- 
tems are set by their cold gas accretion rates (e.g jDave et alj 
l2010t fBouche et al. 2010). In addition, since the average neu- 
tral hydrogen density in the universe has dropped by only 50% 
over 10 billion years in the face of much higher gas consump- 
tion rates, the co ld gas supply must be get t ing continuously re- 
plenished (e.g.,[B auermei ster et al.ll2010t iProchaska & Wolfel 
2009i lPutman et alj|2009al) . 

Despite its importance, the mechanism by which gas is 
replenished in nearby galaxies is a ma jor puzzle (see, e.g., 
ISancisi et af]l2008HPutman et al.ll200 9a. for reviews). Clouds 
of neutral hydrogen are observed around the Milky Way, M3 1 
and a number of other nearby galaxies (e.g..lBlitz et alJll99S " 
IWakker et ail 120071: iMliler et all 120091 iThilker et alj 12004 
but both the fate of these clouds and the fraction of their 
gas that they deposit in their host's disks are uncertain. Di- 
rect searches in nearby groups of galaxies have failed to 
uncover massive po pulations of clouds of neutral hydrogen 
dPisano et alj 120071) . Nevertheless, estimates for the accre- 
tion rate due to high velocity HI clouds (HVCs) have been 
made for a handful of galaxies, and typical values range from 
10% to 20% of th e current SFR of the parent galaxy (see 
Sancisi et al. 2008, for a review, and table|2]below). 

The accretion of gas rich satellite galaxies provides a sec- 
ond potential sour ce of fresh gas for massive star forming 
host galaxies (e.g.. ISancisi et al.ll2008l) . Larger satellites have 
the potential to carry fresh gas to the star forming disk of 
their hosts when they merge, and such a scenario is likely 



in the case of M31 and M33 for example (Putm an et alJ 
l2009d) . However, due to the gas stripping that occurs dur- 
ing the most common minor mergers, satellites will tend to 
deposit their gas at large r adii, well beyond the radius of 
their host's star form i ng dis k (lPeekl2009l : lGrcevich & Putmanl 
120091) . ISancisi et alJ (120081) estimate that some nearby galax- 
ies accrete gas from satellites at a rate of ~ 10 - 20% of 
the SFR of their hosts, alt hough this number cou ld be sig- 
nificantly lower on average (Kauffman n et alj|2010l) . 

Such low gas accretion rates could imply that galaxies are 
presently in the stage of exhausting their gas reservoirs and 
ceasing star formation activity. However, this would leave 
questions about the gas supply during the past several billion 
years, because the evolution of the disk galaxy population is 
inconsistent with a comp l ete shutdown of star formation (e.g., 
iBauermeister et alj|20iq iBell et alj|2007l) . 

So far, predictions of theoretical models as to how gas 
is delivered to galaxy disks are quite uncertain. While 
theoretical studies predict broadly that the average accre- 
tion rate of baryonic mass onto dark mat ter halos should 
be several times their SFRs today (e.g., iDekel et ail 120091 
eq 1), the timescale and route by which the baryonic mass 
reaches galaxy disks thems elves is difficult to predict (e.g., 
iFaucher-Giguere et alJ201 A) . Most of the baryons accreted at 
low redshifts are thought to encounter previously accreted ma- 
terial and suffer shocks, i nstabilities, and fragmentation (e.g., 
iKeres & Hernquistl20 09): these events may mix the accreting 
gas with hot halo gas before it has a chance to reach the star 
forming disk where it is needed. Although such hot gas may 
cool at the disk-halo interface and contribute to star formation 
at a later time, this gas may not be immediately available. 

Given that a large gas reservoir may surround galaxy disks, 
an inviting solution to the gas deficit problem is to tap this 
hot gas through mixing with colder gas on hot phase/cold 
phase boundaries. This mixing would elevate the cooling 
rate of the hot coronal gas, and allow it to sink onto the 
galaxy, thereby producing net accretion . Such mixing could 
occur in the galactic fou ntain (e.g., Mari nacci et alJ 120101 : 
iFraternali & Binnevl 120081) and, intriguingly, might explain 
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the decrease in angular momentum that is observed in gas 
with increasing height above the Milky Way disk. Mixing has 
also been proposed to cool gas out of the ha lo at the interface 
betwe en the hot halo and dwarf satellites (Bland-Hawthorn 
2009). Unfortunately, the predictions of these models are 
rather uncertain because magnetic fields, conduction and 
poorly constrained gas density, and temperature, can all have 
leading-orde r effects on the final mass of gas that cools durin g 
mixing (e.g jMarinacci et al.ll2010tlHeitsch & Putmanl2009l) . 
Even without physical uncertainties, the resolution required to 
track the detailed cooling of the halo gas may be beyond what 
is achievable by current simulations (e.g., an ability to resolve 
the cooling layer may be necessary), so the thermodynamics 
and fate of the accreting gas are both difficult to model and 
are poorly constrained. 

Before exploring unconstrained and difficult to model 
mechanisms for accretion, all other channels should be fully 
understood. A significant internal gas channel - the stel- 
lar population itself - has not been rigorously explored in 
this context. Stellar populations can shed a significant frac- 
tion (~ 30 - 50%) of their mass over the Hubble time 
through secular processes, particularly via shedding of the 
stellar envel opes during the Asymptotic Giant Branch (AGB) 
phase (e.g., [Va ssiriadis & WoodlTT99H iHurley et all l200a 
iGroenewegen et alj|2007f) . In comparison, the Milky Way re- 
quires only ~ 2% of its stellar mass to be returned to ISM per 
billion years to balance its gas consumption by star formation. 
Stellar mass loss should be an accessible fuel for star forma- 
tion in disk galaxies because winds are typically injected di- 
rectly into the relatively dense interstellar medium (ISM) with 
cool temperatures < 3 x 1Q 3 K, l ow velocities ~ 10 km s" 1 and 
largely in molecu lar form (e.g. lKnapp|[l990t lMarengoll2009t 
Libert et al. 2010, and references therein). 

Mass return by stellar populations is noted in the seminal 
study of dRobertsI 1 1 963h . which was the first to discuss the 
transi ence of galaxy p opu lations as a resu l t of ga s consump- 
tion. ISand age ( 1986) and iKennicutt et alj {1994) both made 
the important observation that recycled material could extend 
gas consumption timesca les i n the nearby disk galaxy pop- 
ulation. More recently, iBlitzl (119971) has pointed to observa- 
tions demonstrating that the mass of gas lost in the Milky Way 
could be a significant source of fuel. Nevertheless, in recent 
literature, stellar mass loss is generally not considered as part 
of the gas supply budget outside of some chemical enrichment 
studies. 

Many codes for simulating galaxy formation have only 
recen t ly begun to incorporate stellar mass loss (|Katz et alJ 
Il996t iKravtsov & Gnedinl 12005b IStinson et alJl2006l) and re- 
cent studies indicate that it can be an important physi- 
cal ingredient, significantly influencing galax y morphologies 
dAgertz et al JI20 1 Ot iMartig & Bournaudlr2010 ) and providing 
substantial fuel sou rce for continuing star formation at z — 
(ISchave et alJl2010h . 

The goal of this study is to carefully assess the potential 
importance of stellar mass loss in fueling star formation of 
late-type galaxies at low redshifts (z < 1), by making use of 
observational constraints from star formation measurements. 
To this end, we estimate the total loss rates averaged over en- 
tire galaxies for empirically motivated star formation histories 
(SFHs) and show that stellar mass loss could be an important 
source of fresh gas for star formation at late epochs. In par- 
ticular, we explicitly demonstrate that gas returned to the ISM 
by stars can provide most of the gas required to maintain cur- 
rent level of star formation in a number of nearby galaxies, for 



which the observed accretion rate of halo gas clouds appears 
to be insufficient to fully resupply the gas reservoir. 

To estimate the global stellar mass loss rate for a galaxy 
we need a mass loss model that describes the mass loss rate 
of a single age stellar population as a function of time and 
a star formation history describing the age distribution of all 
stellar populations in a galaxy. Mass loss modeling and its 
uncertainties are discussed in §|2] We calculate star formation 
histories in §[3]based on empirical measurements of the slope 
and evolution of the relation between stellar mass and the star 
formation rate of star forming disk galaxies. The resulting 
global stellar mass loss rates that our model predicts are pre- 
sented in §[U where we also derive mass loss rates for several 
nearby galaxies that are then compared to the observed differ- 
ence between star formation and cold gas infall. Our results 
and conclusions are discussed in § [5] Finally, the Appendix 
explores the extent to which recycled material might accu- 
mulate over time based on the results of cosmological galaxy 
formation simulations and simple models that reproduce those 
simulations. 

2. MODELING MASS LOSS 

To first order, the mass recycled by a single-age stellar pop- 
ulation (SSP) comes from mass that is lost by stars that have 
evolved through the luminous, wind-driving, SNe, red super- 
giant and/or asymptotic giant branch stages, minus the to- 
tal mass in compact stellar remnants. The fraction of mass 
lost by a SSP population at a given time t after its birth is 
therefore approximately the fraction of the population's stellar 
mass that was initally possessed by stars with masses above 
the main sequence turn off m > m to , minus the fraction of 
mass that rema ins locked in the r emnants from those same 
stars (see e.g.. Kalirai et al.l 120081 for empirical constraints 
on the remnant mass fraction). Here we will use mass loss 
rates tabulated from stellar evolution tracks to estimate the 
amount of mass lost by an SSP, b ut calculating mass loss 
from a main sequence turnoff time (Raiteri et al ] 1 19961) and 
the initial-final ma ss relation (from Kalirai et al. 2008) as in 
Agertz et af| (120101) . results in only small differences (5 - 10% 
less mass loss in the first 2 Gyr) from what is presented below. 

The underlying IMF sets the fraction of stars with mass 
above m lo (t). A number of commonly used IMFs are param- 
eterized in table [T] where we define £ to be the number of 
stars per logarithmic mass interval. These parameterizations 
include broken power laws for different mass intervals, 

£(ln(m)) oc m~ r 
and lognormal distributions, 



£(ln(m)) oc exp 



-(logio(w) - log 10 (m c )) 2 



Icr 1 



Corresponding cumulative stellar fractions for different IMFs 
are plotted in the left panel of Figure [1] which shows that 
systematic differences between IMFs at 1M , which corre- 
sponds to the main-sequence turnoff mass after the 10 Gyr 
, are « 20%. At larger masses, corresponding to younger 
population ages, the differences can be as large as a factor of 
* 2 - 3. 

The right hand panel of figure Q] shows cumulative stel- 
lar mass loss as a function of time for SSPs with the same 
IMFs and of solar metallici ty. These were calculated in 
FSPS dConrov & GunnEOlOl) using the Padova stellar evolu- 
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tion models 1 (based on Marigo et al. 2008; Marigo & Girardil 
l2007t iGirardi et al J2000I) 

The figure shows that the majority of stellar mass loss oc- 
curs in the first two billion years, but there is a persistent slow 
rate of mass loss continuing to late times. The differences in 
mass loss for different IMFs after several billion years can be 
as large as a factor of two. However, for the range of IMFs 
considered to be likely for stellar populations in normal galax- 
ies (see, e.g., [Kroupa 2007) the differences are < 30%. 

The fraction of mass lost for ea c h case is well fit by the 
functional form of Jung wiert et al.l (120011) . Namely, the cu- 
mulative fraction of mass lost by a stellar particle at a time t 
since its birth is given by 



/ml(0 



C ln(l + l) 



(1) 



fit parameters Co and A are given in table Q] 

Although the rates of mass loss in any particular phase of 
stellar evolution may be metal dependent and mass loss rate 
measurements for individ ual stars in that phase ma y exhibit 
factor of < 4 scatter (e.g.. iMauron & Jo sselin 2010), the im- 
portant quantity here is the observed initial-final mass rela- 
tion. This appears to be a weak function of metallicity, such 
that aZ = O.1Z population will lock approximately the same 
amount o f mass in remnants as a solar meta llicity popula- 
tion (e.g.. lKalirai et al.ll2008llMarigo & Girardill2007l) . Uncer- 
tainty in the IMF is therefore the dominant source of uncer- 
tainty in estimating the amount of mass loss from a SSP. 

To re present the plausible range of IMFs (see recent re- 
view by [Bastian et al] 120101 for an extensive discussion of 
observational constraints on the possible range of IMFs), we 
adopt the commonly used Chabrier ( 2003, hereafter C03) IMF 
and its steep version, stpC03. The latter has the same log- 
normal form a s the C03 IMF below 1M Q , but a slope that 
is steeper than Salp eterl (119551) by 0.35 above 1M . These 
choices bracket the range of IMFs that are favored by ob- 
servations 2 : they include the observed flattening of the mass 
funct ion at low masses ( iMiller & Scalol [l~979t iRroupa et al.l 
1993), and also encompass the range of high mass slopes that 
are allowed by observatio ns, when theoretical uncertainties i n 
unresolved binary stars (Kroupa et "aD 1 1 99H IRroupa! 120071) . 
dynamic al evolution in clusters and other effects are con- 
sidered dKroup][&Weidni] |200l iBaitiin^tai] |20i3) . Al- 
though they are listed in our tables and shown for compari- 
son, we choose not to brac ket the low mass loss case w ith ei- 
ther th e original power law Salpeterj (119551) IMF or the lScalol 
( 1986) IMF becaus e both are thought to predict too many low 
mass stars (< 1 M [M iller & Scalolll979t iKroupa et aDll993b 
Scalo 1998; Chabrier 2001), and over predict dynamical mass 
dBell & de Jongj|2001[; ICappellari et al.ll2006l) . Unless other- 
wise noted, C03 will be used as the fiducial IMF in this study. 

3. THE STAR FORMATION AND MASS LOSS HISTORIES OF 
GALAXIES 

With mass loss rates for SSPs in hand, the age distribution 
of stellar populations that are losing mass, i.e. star formation 
histories, are needed to determine the relevance of the galaxy- 
wide mass loss to the gas budget and for fueling continuing 
star formation at low z. 

1 http: //stev.oapd. inaf . it/cgi-bin/cmd 

2 The C03 and stpC03 cumulative mass functions vary by no more than 5% 
from the Kroupa 12001) and Kroupa et al. 1 1993) mass functions respectively, 
meaning their mass loss rates are very similar. 



To this end, we will use a model of SFHs derived from the 
observed relation between SFR and stellar masses of galaxies 
at different redshifts. The SFHs can be derived in this way 
because stellar mass in disk galaxies is primarily the integral 
of the in-situ SFR modulo mass loss. Given the observed evo- 
lution of the SFR-M, relation with redshift for actively star 
forming galaxies, SFHs for galaxies of a given current stel- 
lar mass, M*o, can be derived by starting with M,o, integrat- 
ing the stellar mass back in time and determining the corre- 
sponding SFR at each epoch z using the evolving M*(z) and 
the SFR-M, relation at that epoch. Such a procedure is justi- 
fied because scatter in the SFR-M, relation is relatively small: 
the distribution of SFR(M,) at a giv en g is approximate ly log- 
normal with a lcr, of a 0.3 dex dNoeske et alj|2007l) . This 
implies that the galaxy population spends « 95% of its time, 
and builds up « 87% of its stellar mass, during periods when 
SFRs are within a factor of four (< 2cr„) of the median at a 
given epoch. Nevertheless, we explore the effect of the scatter 
i n the SFR-M, relation on mass loss rates more thoroughly in 
§ 13.31 A more detailed discussion of average star formation 
histories derived from measures of instantaneo us star forma- 
tion rate will be presented in a separate paper dLeitnerll201 ll 
in preparation). Although SFHs of individual galaxies might 
be determined more accurately by fitting population synthesis 
models to galaxy spectra, the approach we adopt here should 
provide an average estimate of the importance of stellar mass 
loss for galaxies of a given stellar mass. 

Implicit in this strategy for deriving SFHs is the assumption 
that mergers play no significant role in stellar build-up and 
that a galaxy's stellar mass changes solely because it forms 
stars or loses mass through stellar mass loss. Mergers have 
the potential to split M, as lookback time increases, thereby 
slightly increasing the SFR for the ensemble of progenitors 
compared to the values derived using our model. If mergers 
were an important factor in the growth of low-z galaxy stel- 
lar mass, the SFR derived using our model would be under- 
estimated resulting in galaxies forming more of their stellar 
mass at shorter lookback time. There is evidence, however, 
that the effect of mergers on the buildup of stellar mass is 
not large for most disk gala xies of M,q $ 1 n Mp) at low 
reds hifts. For example, both iBell et all (120071 see Figure 4) 
and iDrory & Alvarezl d2008l) show that a scheme similar to 
ours reproduces the buildup of the stellar mass function of all 
galaxies w ith only a minor contributi on from mergers. Fur- 
thermore, IConrov & Wechslerl (120091) derive average SFHs 
for galaxies using the subhalo abundance matching method 
and show that mergers of galaxy halos cannot play a signifi- 
cant role in the evolution of the stellar mass function of galax- 
ies at z < 2. 

3.1. The Observed Evolution of the SFR-M* Relation 

A number of studies have explored the evolution and 
shape of the SFR-M , relation at complementary UV and 
infrared wavelengths dNoeske et al.ll2007b lElbaz et al] 1 2007 1: 
I Salim et al l l2007t iDaddi et al.ll2007l l2009t IStark et alj|2009 



Gonzalez et al.l 20101 see t he compilation of results by 



Dutto netal.H2010h . Recently. lOliver et al.l d2010l) performed 
a far-infrared stacking survey at 70 and 160 pm, and included 
100 um and 160A»n obs ervations from the Herschel satellite 
(Rodi ghiero et al.l 120101) . which allowed them to probe re- 
gions of the SEP l ess aff ected by uncertainties in dust absorp- 
tion. lOliver et"ai1 d2010l) derived a fit to the median evolution 
and shape of the SFR-M, relation of star forming galaxies 
in their sample. The redshift evolution they derive is plot- 
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Figure 1. Left panel: The fraction of stellar mass in stars of mass > m for a number of commonly used IMFs. Right panel: The fraction of mass lost as a 
function of age for a single-age stellar population. Chabrier()3 and Chabrier03 Steep were chosen to bracket the effect of likely IMF choices on mass loss rates in 
our study (see text for details). 



Table 1 

Mass loss and the IMF 



IMF 


Functional form 


Mas 
Co 


s loss fits 
A[yr] 


Chabrier (20031 TC031 
(similar to Kroupa 2001) 


m < 1M Q : log-normal (<x = 0.69, m c = 0.08) 
m > 1M Q : power-law (T = 1.3) 


0.046 


2.16 X 10 5 


Steep Chabrier f2003) [stpC03] 
(similar to Kroupa et al. 1993) 


m < 1M Q : log-normal (<x = 0.69, m c = 0.08) 
m > 1M Q : power-law (r = 1.7) 


0.051 


1.33 x 10 7 


Scalo (1998) TS981 


m < lM e : H = 0.2 
1M <?n< l0M e : T 2 = 1.7 
m > 10M G : T 3 = 1.3 


0.062 


6.97 x 10 6 


Miller & Scalo (1979) 


log-normal (<t = 0.74, m c = 0.95) 


0.058 


6.04 x 10 6 


Salpeter(1955) 


T = 1.35 


0.032 


5.13 x 10 5 


Scalo (1986)" 


m < 0.2 : Ti = -2.6; 0.2 < m < \M : V 2 = 0.80 
IM B <m < 10M o : T 3 = 2.25 
m > 10M o : T 4 = 1.45 


0.055 


1.25 x 10 8 



"A fit from Mo et al. 12010) with an extension to m < 0.2. 



ted as a solid line along with the points from previous stud- 
ies in Figure [2] The figure shows a general consistency be- 
tween different observations at z < 1 (points have been scaled 
to M» = 1O 1O - 75 M according to the SFR-M, relati on pro- 
vided for each separate data set in Dutt on et alJu OlO). Stud- 
ies including quiescent galaxies (e.g., Zheng et al. 2007|) also 
match the recent IR-derived results dOliver et al.| [2010). indi- 
cating that systematic errors are under control (although see 
Rodig hiero et al.l 120101 for recent L^g-S FR calibration from 
Herschel). In this study, we use the fit of lOliver et ail (l20loh 



at z < 2. 

At z > 1 uncertainty in the normalization and slope 
of the SFR-M* relation increases and there is a possibility 
that an evolving IM F may affect the evolution dDavel 120081 : 
IWiMns et alll2008l but see lBastian et al.ll2010l) . Fortunately, 
the build-up of stellar mass at z > 2 does not significantly af- 
fect estimates of stellar mass loss at z = for late type galax- 
ies, which are the focus of our study. Furthermore, late-type 
galaxies with M, < 1O U M appear to build up most of their 
stellar mass at z < 1, so the evolution of star formation rates at 
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Figure 2. Observations of the evolution of the normalization of the SFR-M. 
relation at the stellar mass M, = 10 10 - 75 M R w ith redshift for various data 
sets. 



The red points are fr o m the IDutton et alj 120101) comp i lation of dat a 
from IGonzalez et alj f20TOh; IStark et alj 12009); Daddi et al. 1 20091 120071) ; 
lElbazetalJ 12007): Noes ke et alj <200 7): Salim et alj <2007D The error bars 
show the estimated factor of two scatter of galaxies about if/(M„,z)- At z £ 2, 
the blue line is the fit to the far infrared data of late-type star forming galaxies 
from Oliver et al. 1 2010), which shows good agreement with previous studies. 
At z > 2, we assume that normalization of the SFR-M, is constant (dashed 
line). 



high redshifts is not very important (see figure[3]l. We mimic 
the high redshift evolution with a redshift independent SFR 
that forms a broken power law with Zbreak = 2. 

Recent measurements of the slope [} in the relation 
SFR oc Ml +l at fixed epoch, show -0.5 < < at 
Z < 1 . Whether /? remains constant with redshift ( Zheng et al. | 



20071) . evolves to a shallo wer value (e.g.. iBrinch mann et al. 
20o4 iNoeske et all I2007t Elbaz et al l 120071). or steep ens 
hiero et al.l2010tlErb et al.l2006tlohver et al. 



(Rodigh 



,2010), 



is 

unclear at present. Here we assume a constant model with 
yS = -0.25. The resulting global scaling we use for z < Zbreak 
is, 



<A(M t ,z)=A (z+ir 



M, 



10 10 - 75 M o 



(2) 



where we adopt the power law slope of a « 3.4 for late-type 
galaxies of M, = 1O IO75 M dOliver et alj|20!ob . We label 
this fit to the median instantaneous SFR measurements i// to 
distinguish it from the SFR that a specific galaxy experiences 
as a function of redshift (its SFH). The normalization Aq to 
(M„ z) (plotted in figure|2]for M, = 10 u M o ) is 2.8M yr-' 
Oliver et al. 201 0). U ncertainty in a and f3 will be discussed 
further in section [3~2l 

For Zhreak < z < 6 we use a — 0, and assume that all initial 
stars in a galaxy form in a burst at z = 6 (these stars constitute 
only a few percent of M,o). 

3.2. Calculating Star Formation and Mass Loss Histories 



(C 



The overall global mass loss rate (GMLR) of a galaxy is 
equal to its SFR convolved with the fractional mass loss rate 
f m i (whose integral is given by eq.Q}: 



GMLR(f) 



Jt 



SFR(f')./ml(f - t')dt' 



(3) 



In our toy model, the evolution of the stellar mass of a 
galaxy is given by, 



M,(f) = i/<M„z) - GMLR(f) 



(4) 



with boundaty condition M,(f = 0) = M»o, where tfr(M t ,t) is 
the instantaneous SFR function from eq.|2] which is a function 
of the evolving mass of the galaxy M». The SFH of the galaxy 
is then given by, 



SFR(f) = ij/(MJt),z) 



(5) 



The GMLR(f) term introduces a complication in that it is 
defined as a convolution of the SFR of a galaxy with its 
fractional mass loss rate up to time t, but the SFH of the 
galaxy is not, at first, known. As a first guess, we can take 
GMLR(f) = 0, derive M„(f), plug the resulting SFH into eq.|2] 
and iterate to convergence, 

M?\t) = iftM?\t) 

AO)/ 



SFR (0) (f) = ip(MT{t),z) 



M { - \t) = i//(M ( t l \ t) - GMLR(SFH (0) , f) 



(6) 



Star formation and stellar mass loss histories are plotted in 
Figure|3]for several representative values o f M*p and the fixed 
median SFR scaling relation discussed in § 13.11 

In the right panel of Figure[3] we show the effect of the un- 
certainty in the parameters a (o~ a ~ +0.5) and ft {a-p a +0.25) 
on the SFH and t he mass loss rates (based on estimates by 
lOliver et"aT1 120 1 01 and dispersion in the literature discussed 
above). The observational lcr uncertainty in the slope and 
evolution of the SFR-M, relation only translates into uncer- 
tainty in the mass loss rate of < 5% at z = 0, although the 
differences can be relevant to the stellar mass loss history. 

3.3. The Impact of Scatter in the SFR-M* Relation 

Although the observed scatter in the SFR-M, relation is 
relatively small, deviations from the median relation can be 
significant for some galaxies and may thus affect conclusions 
about the importance of stellar mass loss. In this section, we 
will therefore consider the effects of populating the scatter in 
SFR-M, relation, with a range of possible SFH scenarios. We 
parameterize the scatter around the star formation rate func- 
tion 4>(M,,z) that an individual galaxy experiences, with the 
time scale, ta, over which deviations around the median rela- 
tion occur, and by the magnitude of scatter, cr,. We assume 
that the scatter is described by a log-normal function with the 
median described by the relation in equation|2] We thus mul- 
tiply iff(M*,z) at a given redshift by by N(t, Ta, 0"*) - a log- 
normal random number sampled every ta for t > t&. To force 
the model to trend toward the assigned SFR(z = 0), we set 
N(t, t a , cr,) = SFR(z = 0)/A for lookback times < t a . 

The effect of scatter is illustrated in Figure [4] which shows 
relative global stellar mass loss rates as a function of assumed 



6 




Figure 3. Left Panel: SFHs (thick lines) for galaxies of different z = stellar mas s, M, : 10 1I M o (red solid line), 10 105 iW o (green dashed line), and 10 100 M o 
(blue dotted line). SFHs here are derived from the observed tf/(M,,z) (see jj |3.2l for details). The thin lines of the corresponding color and line type show the 
global stellar mass loss rates for these SFHs. Right Panel: SFHs (solid lines) and stellar mass loss rates (dashed lines) calculated for a galaxy of z = stellar 
mass of M„o = 10 10 75 M Q with the values of parameters of if/(M,,z), a and f} (see eq.ff} varied by lcr (cr a ss 0.5 and crp rj 0.25) from their average measured 
values (lines of different color, with color correspondence indicated in the legend). Note that stellar mass loss rate at low z is insensitive to lcr variation of the 
parameters of ip(M„ , z) at this level. 



ta for model galaxies with M* = 10 ia75 M G (the stellar mass 
of the Milky Way) and SFR(z = 0) corresponding to the 
median SFR at that mass, as well as star formation rates at 
+ lcr t from the median. Mass loss rates are normalized by 
the population-average rate in the ta = case. The overall 
amplitude of each line in the vertical direction therefore in- 
dicates the factor by which the mass loss rate is biased with 
respect to the case wherein galaxies, at any instant, draw their 
SFR randomly from the lognormal population-wide probabil- 
ity distribution function of SFRs. The shaded regions covers 
the area between the 16% and 84% outliers from the median 
mass loss rate at a given {SFR(z = 0), ta}. 

The shaded region in the figure shows that that scatter in the 
SFR-M* relation results in small scatter in the mass loss rate 
for a given galaxy. The different present day SFRs also result 
in an overall change of the mass loss rate, although this effect 
disappears for very short t&. 

Galaxy-wide stochasticity would presumably have a cata- 
lyst, such as a merger, tidal interaction, or rapid increase in 
gas accretion rate. These processes should operate on the 
crossing time of the halo (~ 1 Gyr) so we adopt ta = 500 Myr 
as a fiducial minimum value. We do not enforce a maximum 
ta, thus we use ta = 14 Gyr for the scenario in which the 
scatter around the median relation for the disk population is 
not generated by stochasticity in individual galaxies, but is in- 
stead caused by persistent environmental effects; these effects 
would be different for different galaxies in the population, but 
similar over time for a given galaxy. In other words, in such a 
scenario, a galaxy that has a SFR below if/(M t ) at some high 



redshift will maintain a low SFR with respect to if/(M*,z) at 
all z. The overall uncertainty in the global mass loss rates for 
a given mass loss model, encompassing both model depen- 
dence and the +lcr A outliers from stochastic variation with 
ta = 500 Myr, can then be quantified as the GMLRs encom- 
passing these two regions. This uncertainty is marked by the 
black vertical dotted lines for each S FR(z = 0) in Figure |U 
and this method will be used to describe model uncertainty in 
the results below. 

The star formation and stellar mass loss histories for model 
galaxies with ta = 14 Gyr are plotted in the right panel of Fig- 
ure!?] for me same representative range of present day SFRs. 
The lower SFR(z = 0) case results in less star formation at late 
times, and hence more star formation at high redshifts, as the 
total stellar mass is kept constant. This push-back also results 
in less stellar mass loss at low z, since a population loses more 
of its mass when it is younger. However, the decrease in mass 
loss rate is not proportional to the decrease in SFR(z = 0), 
as is assumed by the instantaneous recycling approximation. 
This is because old populations continue to shed mass. In fact, 
the mass loss and star formation lines are very close to each 
other for the low SFR(z = 0) case in FigureH] indicating that 
stellar mass loss could be responsible for fueling almost all of 
the star formation in such a galaxy. Conversely, in galaxies 
with a higher current SFR for a given stellar mass loss pro- 
vides a smaller fraction of the fuel required to maintain the 
star formation at the same level. 

4. IMPLICATIONS OF STELLAR MASS LOSS FOR THE GAS SUPPLY 
IN GALAXIES 
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Figure 4. Left Panel: The z = stellar mass loss rate for galaxies of M,q = 10 10 75 M Q as a function of ta normalized to the population's average mass loss 
rate in a scenario where galaxies draw their SFR randomly from a lognormal distribution about i//(M,,z) at all times (i.e., r A = 0, see text for details). Solid, 
short dashed, and long dashed lines correspond to galaxies with present day SFR that are at the median (Aq), 16th (0.5Ao) and 84th (2Aq) percentile in the scatter 
about ijj(M,,z) at M»o = 10 10 - 75 M G . Green regions encompass the 1o"a scatter for a population about their corresponding relative mass loss rate, while blue 
lines show the ta = 14 Gyr limit for each present day SFR. The ta = 14 Gyr cases are systematically lower for each present day star SFR because the mean of 
a lognormal distribution is offset from the median. We take r A > 500 Myr so the difference between the mass loss rate with ta = 500 Myr and 14 Gyr, noted 
with vertical black dotted lines, marks the uncertainty in the mass loss model. Right Panel: Star formation (thick red) and stellar mass loss (thin blue) histories 
atM«o = 10 1075 M Q forr A = 14 Gyr. Line styles are the same as on the left. The model shows that stellar mass loss can provide a larger fraction of the fuel 
required to maintain star formation in galaxies with low SFR(z = 0). In all cases, the global stellar mass loss rate is a significant fraction of the present day SFR. 



4.1. z — Stellar Mass Loss 

Present day stellar mass loss rates for galaxies, calculated 
assuming both C03 and stpC03 IMFs (see § |2]i are plotted 
in Figure [5] as a function of galaxy stellar mass. The addi- 
tional infall that is required to match the consumption of gas 
by star formation is the difference between the SFR and stellar 
mass loss; it is shown by the hatched bands. The width of the 
shaded regions corresponds to uncertainty about the C03 and 
stpC03 IMFs caused by stochasticity (with r A > 500 Myr) 
that can co mpri se the scatter in the SFR-M* relation, as dis- 
cussed in § 13.31 Note that this uncertainty is relatively small 
and does not affect our conclusions about the importance of 
stellar mass loss. 

The middle panel of Figure [5] shows that stellar mass loss 
dominates over infall across the entire plotted stellar mass 
range. For a stpC03 IMF, stellar mass loss still makes up a 
significant portion of the star formation but, on average, infall 
should comprise 45 - 60% of the star formation rate. 

The side panels of Figure|5]show the amount of stellar mass 
loss for galaxies with present day star formation rates +lcr t 
(i.e., a factor of two) away from the median SFR-M* relation. 
The left panel shows that for galaxies with SFR below the me- 
dian, stellar mass loss dominates over infall for both C03 and 
stpC03 IMFs. The right panel shows that galaxies with SFR 
higher than average require more significant infall to maintain 
their SFR. For the stpC03 IMF infall should contribute about 



60% of the gas to maintain current SFR, whereas for the C03 
IMF, stellar mass loss still dominates, but infall must make up 
for about 40 - 50% of the gas required to maintain the SFR. 

4.2. Recycling Epoch 

The epoch at which the stellar mass loss rate becomes larger 
than half of star formation rate is plotted in Figure[6] This pe- 
riod can be thought of as the recycling epoch, when stellar 
mass loss provides most of the new gas needed to replenish 
what is consumed by star formation. The figure shows that 
the recycling epoch occurs at z ~ 0.7 - 0.9, z ~ 0.3 - 0.6 and 
~ 0.05 - 0.2 for the galaxies with respectively systematically 
low, median, and high star formation rates relative to i[f(M*,z)- 
The recycling epoch is sensitive to the assumed IMF; for ex- 
ample, assuming stpC03 IMF instead of C03, galaxies that 
follow n/(M*,z) enter the recycling epoch at z ~ 0.1. 

4.3. The Gas and Star Formation Budget of Nearby Galaxies 

It is only possible to estimate the cold gas accretion rate 
for a handful of galaxies. Nevertheless, available estimates 
are all considerably smaller than the current star formation 
rates of host galaxies, a fact that has been used to raise 
concerns of a siz a ble gas deficit jFraternalil 120 Id l2009t 
I Sancisi et alj|2008t iFraternali & Binnevl 120081 iPutman et all 
2009b; Marina cci et alj|2010l) . It is thus interesting to evalu- 
ate the potential contribution of stellar mass loss. 
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Figure 5. Average star formation rate (red dashed lines) and global stellar mass loss rate (blue dashed lines) as a function of stellar mass. The green solid lines 
show the difference between SFR and stellar mass loss rate. Lines without points are calculated assuming no stochastic scatter (i.e., = 14Gyr), and lines with 
points are for the models with scatter assuming = 500 Myr timescale with a C03 (triangles), and stpC03 (squares) IMFs. Shaded regions represent modeling 
uncertainty for the C03 (lightly shaded) and stpC03 (densely shaded) IMFs, which includes plausible range of stochasticity in the SFR-M, relation. The central 
panel is for median present day star formers, while the left and right panels are for galaxies with star formation rates higher and lower than the median by 1<t», 
respectively. Note that in all panels stellar mass loss rate is a significant fraction of the SFR and is comparable to the SFR in the left panel (i.e., low SFR case). 



12 



10 



O 
CO 
Xl 

o 
o 



8 - 



2 - 



T 



SFR(z = 0) = 0.5f(M.,z) 

.SFR(z = 0)=f (M.,z) 
.SFR(z = 0) = 2V(M„,z) 




_l I I I I I I L 



z = 3 
2 



10 10 



5xl0 10 10 11 



M. [MJ 



Figure 6. The lookback time at which the global stellar mass loss rates be- 
come larger than half the star formation rate, assuming a C03 IMF. The short 
(long) dashed line show such time for galaxies in the 16th (84th) percentile 
in SFR(z = 0), while the solid lines show results for the median SFR-M» 
relation. 



The gas budgets for the galaxies with available estimates 
of gas infall rates are presented in Table [2] and Figure [7] in- 
cluding the expected contribution from stellar mass loss. To 
estimate the mass loss rate we use observational constraints 
on SFR(z = 0) and M*o as input and calculate present day 
stellar mass loss rates through modeled SFHs, as described 
above. For galaxies with a range of SFR(z = 0) quoted in 



the table, we use the high value to calculate the mass loss 
rates. This choice is conservative, because mass loss will 
make up a smaller fraction of the gas supply in galaxies with 
higher SFR at fixed stellar mass, as was shown above (see, 
e.g., Fig.@]right panel). Ranges quoted for global stellar mass 
loss rates incorporate t he s catter and modelling uncertainty 
that was quantified in § 13.31 

The table and figure show that total gas supply rates (stellar 
mass loss + gas infall) are close to star formation rates. Con- 
servatively, the factor of several discr epancy between star for- 
mation and cold-HI gas supply (e.g., Fraternali 2009) should 
reduce to a factor of two discrepancy for some galaxies, and 
to zero or a surplus in others. 

Furthermore, a less conservative accounting of infall may 
include factor of two or larger corrections to infall rates from 
(1) a 33% correction for the helium fraction, and (2) the 
fact that neutral gas only accounts for a fraction of the in- 
falling cloud mass. The neutral fraction in the HVC clouds 
is not very well constrained at present with observational es- 
timates ranging from ~ 10 - 20% for HVCs that are distant 
dMalonev & Putmanll2003b to ~ 50% for the clouds near the 
disk (IWakker et alj|2007t iHill et al.ll2009l) . Despite the uncer- 
tainties, it is clear that ionized ga s constitutes a non-ne gligible 
mass fraction of the halo clouds (ISembach et al.1 12003) and so 
the accretion rate observed in HI can only be a lower limit. 
There could also be additional possible sources of fresh gas, 
such as infall from gas rich satellites (0.1 - 0.2M o yr _1 ) and 
accretion tha t is lost in con f usion noise (esp. forNGC 891), as 
discussed by lSancisi et al.l (120081) . Observations of the Milky 
Way suggest that Q.lM^yx " 1 of additional infa ll may have 
been already been observed (Wakker et al. 2008). 

Finally, we note that empirical estimates of stellar mass loss 
exist for some nearby dwarf galaxies, such as LMC, WLM, 
and IC 1613, in which complete samples of AGB stars can 
be constructed and mass loss for in dividual AGB stars d e- 
termined from the IR observations (iJackson et al.l l2007alfbl : 
iMatsuura et"aTll2009t ISrinivasan et all 120091) . Observational 
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Figure 7. The gas budget for each galaxy labeled on the bottom axis. Bars on the left of each galaxy denote the observed gas consumption from star formation 
(dark red and magenta). The two sets of bars on the right denote gas supplied from observed infall of new gas (green) and stellar mass loss (blue). The two gas 
supply rate columns are labeled, below, by the C03 and stpC03 IMFs on which they are based. The lighter shades of each quantity indicate the space between 
low and high estimates of that quantity. For stellar mass loss, these are based on the uncertainty estimates from this study; for star formation and infall, they are 
based on the size of cited error bars and scatter between different measurements. Specific values, references and notes are given in Tableff] 



Table 2 

Galaxy properties and the gas budget. 















Stellar Mass Loss Rates 




Galaxy 


M, a a 


Vflai 


median(SFRo) 


SFRo" 


Accretion rate* 


C03 


stpC03 


References 




(10 9 M o ) 


(kms" 1 ) 


(M yr-') 


(M yr-') 


(M yr-') 


(Moyr- 1 ) 


(Moyr 1 ) 


[M.0][SFR ][Accr.] 


M33 


4.5 


104 


0.42 


0.55 - 0.7 


0.05 - 0.8 r 


0.34-0.37 


0.22-0.25 


[1][2,3][4] 


NGC 4559 


6.8 


123 


0.58 


0.6 


<0.05 


0.36-0.39 


0.26-0.29 


[5][6][7] 


NGC 2403 


5-14 


134 


1.0 


0.38-1.3 


0.1 


0.72-0.79 


0.50-0.58 


[8,9][10,11][12] 


Milky Way 


55 


220 


2.8 


0.7-1.45 


0.2 


1.26-1.67 


0.99-1.40 


[13][14,15][16] 


NGC 891 . 


100 


230 


4.3 


3.8 


0.2 d 


2.88-3.30 


2.19-2.62 


[17][18][19] 


M31 


100 


260 


4.3 


0.4- 1.0 


0.0 


1.30-2.30 


1.10-2.10 


[20] [2 1,22,23] [24] 


NGC 2997 


160 


226 


6.2 


5.0 ± 1 


1.2 


4.22-4.75 


3.30-3.80 


[25] [26] [27] 


NGC 5746 


160 


320 


6.2 


0.8 ±0.2- 1.2 


0.2 


1.90-3.44 


1.59-3.14 


[28][29,30][31] 



"Systematic and random errors in stellar mass and SFR are likely < 0.3dex iMo et al. 2010), but difference between stellar masses and SFR derived with C03 
and stpC03 should be < 10% (approximately the difference between the Kroupa et al. 1993 and Kroupa 2001 IMFs found by Borch et al. 2006). * Accretion rates 
are probably lower limits to cold-gas accretion, and do not include helium or ionized gas fractions unless noted. c Infall of surrounding gas clouds at 100 km s 
is assumed, the low (high) value excludes (includes) an ionization correction; ''if all HI in the halo of NGC 891 were to fall toward the disk it would supply 

~ SOAfoyr- 1 ; 

References: [l] [Corbelli (2003); [2] Magrini et al. (2007); [3] Blitz & Rosolowskv (2006); [ 41IGrossi et al.1 j2008l): [51 IBarbieri et at] f200"l. from M,/L k ; [6] 
Fraternali (2010); [7] Barbieri et al. (2005); [8] Lerov et al. (2008); [9] de Blok & McGaugh 1 1996); [10] Lerov et al. (2008); [11] Kennicutt et al. (2003); [12] 
Fraternali et al. (2002); [13] Flvnn et al. (2006); [14] Robitaille & Whitney (2010); [15] Murray & Rahman (2010); [16] Wakkeretal. (2007); [17] Rhode et al. 
(2010); [18] Popescuetal. (2004); [19] Oosterloo et al. (2007); [20] Chemin et al. (2009); [21] Walterbos & Braun (1994); [22] Barmby et al. (2006); [23] 
Williams (2003); [24] Chemin et al. (2009); [25] Louise & Marcelin 1 1983); [26] Hess et al. (2009); [27] Hess et al. (2009); [28] Rasmussen et al. (2009); [29] 
Rasmussen et al. (2006); [30] Pedersen et al. (2006); [31] Rand & Benjamin (2008); 



estimates of total mass loss from the AGB stars for these IC 1613) of their total SFR. Given that additional mass loss 
galaxies range from ~ 10% (LMC) to ~ 50 - 100% (WLM, can be expected from supernovae and red giant stars, these in- 
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dependent, direct estimates are consistent with our model cal- 
culations and also suggest that stellar mass loss may provide 
substantial fuel for ongoing star formation in these galaxies. 

5. DISCUSSION AND CONCLUSIONS 

We have used a simple, empirically-motivated model for 
deriving star formation histories of galaxies of different stel- 
lar masses along with mass loss rates from stellar evolution 
code to show that stellar mass loss can be the most important 
source of gas for fueling continuing star formation in late type 
galaxies at low redshifts (z < 0.1 - 0.5). We also explicitly 
show that the gas supplied by stellar mass loss is comparable 
to the perceived deficit of gas required to fuel continuing star 
formation in a number of nearby galaxies. This conclusion is 
consistent with direct observational estimates of stellar mass 
loss contribution of gas to the ISM of nearby dwarf galaxies, 
such as WLM and IC1613. We presented tests of the effects of 
scatter and uncertainties in the SFR-M* relation on our results 
and have show that our conclusions are not affected. 

We believe that our calculations provide a conservative es- 
timate of the expected stellar mass loss rates. First, they rely 
on star formation histories derived from the observed evolu- 
tion of the SFR-M, relation. Thus, only if current observa- 
tional results on the evolution of this relation are significantly 
biased, would our conclusions be biased as well. Moreover, 
in the Appendix we show that any feedback that delays con- 
sumption of material lost by stars will only further enhance 
the rate that recycled gas gets reprocessed at late epochs. 

Maintaining a gas supply for star formation at z ~ pri- 
marily through recycled stellar stellar mass has implications 
for the recent chemical enrichment histories of galaxies. In 
the case of the Milky Way, ongoing infall of pristine gas 
is thought to be required to match the observed enrichmen t 
history and metallicity distribution (e.g.. lSancisi et al.1 2008). 
However, observations of the Milky Way and other nearby 
galaxies show that the metallicities of planetary nebulae are 
generally lower than or equa l to the metallicities of the HII 
regions around young stars (Stang hellini et al.ll2.010h . mean- 
ing that gas returned to the ISM by such nebulae does not 
need to be diluted by infall in order to provide fuel for cur- 
rently forming stars. Furthermore, detailed enrichment mod- 
els are actually not s o restrictive at late epochs. For exam- 
ple, |Chiapj2inj| (2008) points out that although chemical evo- 
lution models show a need for large inflow rates during the 
early stages of the evolution of the Milky Way disk, a present 
day infall rate of only « O.45M yr 1 is used in the fiducial 
Chiap pini et al.l (1200 ll) model that is consistent with enrich- 
ment constraints. 

Deuterium abundances can potentially provide the most 
stringent constraint on the amount of gas that stellar mass loss 
can supply to galaxies like the Milky Way. Deuterium is com- 
pletely destroyed in stars, so that the deuterium abundance in 
the ISM of a galaxy relative to its primordial value con strains 
the fr action of ISM recycled throug h stars (see, e.g., | J aeel[ 
[2009L for a pedagogical review and IProdanovic et al.l 120101 
for a recent analysis). However, it is not clear how restric- 
tive deuterium constraints are at present. Deuterium abun- 
dance exhibits variation of a factor of ~ 3 - 4 even within a 
kilop arsec from the Sun (Sonneborn et al. 20 00l: iLinsky et all 
120061) and longer lines of sight indicate lower abundances 
dWood et al.l 120041) . One possible interpretation is that vari- 
ation o f deuterium abunda nce is due to its depletion onto dust 
grains (ILinsky et al.ll2006T) . in which case the largest observed 
abundances should be interpreted as the true undepleted val- 



ues. Such assumption leads to conclusion that most of the 
present I SM gas in the Milky Way should be unprocesse d 
by stars (IProdanovic & Fieldsll2008t IProdanovic e"tai1l2010l) . 
However, while dust depletion interpretation of observed deu- 
terium abundance variation is plausible, some puzzles do re- 
main. For example, the correlation between deuterium abun- 
dance along a line of sight and reddening in the same direc- 
tion, which would be expected within s uch framework, has 
not been detected (Steig man et al.ll2~007l) . This implies that 
abundance variation may also be due to incomplete ISM mix- 
ing, in which case the constraints on recycled fraction can 
be substantially relaxed. Finally, we note also that in the 
context of this study, the deuterium abundance in the solar 
neighborhood may not be representative for the global con- 
straints for the entire Galaxy. This is because most of molec- 
ular gas and star formation in the Milky Way occurs in a ring 
at ~ 3-6 kpc from the Galactic center and the deuterium 
abundan ce in that region is substan tially lower than around 
the Sun dLubowich & PasachoflfeOlOl) . 

An interesting implication of our results is that stars can be 
forming largely from recycled stellar material for cosmolog- 
ically significant periods of time. The question is then why 
some galaxies with significant stellar mass do not form stars 
during late stages of their evolution. Alth ough some cold 
gas is observed in early type galaxies (e.g.. iBouchard et"al] 
120051; IWelch et"ail 120101) . most elliptical galaxies have star 
formation rates much lower than the expected rate of stel- 
lar mass loss. Although detailed analysis of this issue is be- 
yond the scope of this study, the likely explanation is that gas 
lost by stars in spheroidal systems orbiting with high veloc- 
ities within tenuous hot halo is efficiently mixed with halo 
gas. The situation is an alogous to a well -known "cloud in 
the wind" problem (e.g., Agertz et al. 2007) with hot halo gas 
acting as high-velocity wind blowing past a red giant enve- 
lope. The cloud is expected to be unstable to the Kelvin- 
Helmholz (KH) instability and to be disrupted and mixed 
with hot gas on a short time scale . The extent of this mix- 
ing is a subject of active r esearch (jBregman & P arriotll 120091 : 
iParriott & B regman 2008). At the center of early types, densi- 
ties of recycled material should become so high that only reg- 
ular AGN activity may keep recycled material from condens- 
ing (ICiotti & Ostrikej] 120071) . Interestingly, the abundance 
of dust expected to come from the envelopes of red giant 
stars according to stellar mass loss mo dels is also not ob- 
served in nearby globular clusters (e.g., iBarmbv et al.l 120091; 
iGnedin et al.ll2002t Ivan" Loon et alJl2009h . In this case, globu- 
lar cluster stars are also moving through Milky Way halo gas 
at high velocity, so no accumulation of cold gas shed by stars 
may be possible as that gas is, again, experiencing ram pres - 
sure, KH instabilities, and mixing (e.g. . Ivan Loon et al.i 2009). 

The situation is likely to be different for disk galaxies be- 
cause disk stars are releasing their gas directly into co-rotating 
ISM, rather than rapidly flowing hot halo gas. The lost gas 
can therefore join ISM and become readily available for star 
formation. Therefore, our results and conclusions about im- 
portance of stellar mass loss should be primarily applicable to 
late-type disk galaxies. In this respect, it is interesting that ob- 
servations indicate that there is more gas in the luminous SO 
galaxies than in the elliptical galaxies of similar luminosity 
(IWelch et al.ll20loh . 

Even in the disk systems, however, the stellar mass loss may 
only be efficient in fueling star formation in the presence of 
at least some gaseous disk. Gas lost by stars in a pure stel- 
lar disk (without any corotating gas) would be subject to the 
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same efficient mixing with halo gas, as stars that are moving 
with high velocities through slowly rotating halo. Therefore, 
stellar mass loss may become inefficient in fueling continuing 
star formation when gas surface densities falls to some small 
value. 

The significant role of stellar mass loss in fueling contin- 
uing star formation during late stages of evolution of galax- 
ies implies that it is very important to include this process in 
cosmological simulations of galaxy formation. Inclusion of 
stellar mass loss in simulations is straightforward: for a given 
choice of IMF, the mass of gas returned to ISM in a given in- 
terval of time is given by the differential of equation Q] with 
appropriate parameters (see Table[TJ. This mass can be added 
directly to a host cell of stellar particle in a grid simulation or 
distributed over neighboring gas particles in SPH simulations. 
Simulations with and without stellar mass loss can be used to 
explore other possible e ffects of this process, suc h as effect 
on galaxy morphologies (Martig & Bou rnaudl2010l) and their 
stellar masses. The latter will depend on the star formation 
histories of galaxies and thus may exhibit systematic varia- 
tion with host halo mass and contribute to the scatter of stellar 
mass at a fixed halo mass. 
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APPENDIX 
MODELING GAS REPROCESSING 

Until now we have dealt with mass loss from a stellar population at a particular instant, ignoring that gas returned to ISM may 
not be immediately accessible to form stars and may instead accumulate. Using cosmological galaxy formation simulations, we 
now address issue of how much gas might be available from stellar recycling to be reprocessed to form stars at low redshift, given 
gas reprocessing physics over galaxy's mass loss history. 

Simulations 

Cosmological galaxy formation simul ations of MW-sized systems were carried out using Eulerian, gas dynamics+A^-body 
adaptive refinement tree (ART) code (lKravtsovlll999t iKravtsov et al.ll2002t iRudd et al.1 120081) . The simulations followed the 
evolution of dark matter and baryons in a box of comoving size Lbox = 20/z _I Mpc in a flat ACDM cosmology: Qo = 1 - Sl\ = 
0.258, Q. b = 0.044, h = 0.72, n s = 0.96, cr 8 = 0.80 

A standard "zoom-in" approach was used to generate initial conditions in order to achieve high mass resolution, while adaptive 
mesh refinement was used to increase the spatial resolution in regions where baryons collapse to form galactic disks. Milky Way- 
sized systems were selected for resimulation from a low -resolution non-radiative simulation at z = 0. Particles within three virial 
radii of the around the center of each object were then replaced at the initial epoch (z = 62) with eight times more high-resolution 
particles of mass hjdm = 2.82 x IQ 7 h~ l M Q with a spatial distribution constrained to match the low resolution initial conditions, 
while also including additional smaller scale power. There is no contamination by low resolution dark matter particles within the 
virial radius of our halos at z — 0, but the large-scale tidal field remains accurately sampled using this technique. 

The simulation box was initially resolved by a uniform 128 3 grid and additional successive refinements were introduced in 
collapsing regions that were occupied by the high-resolution dark matter particles. The cells were refined if their gas mass 
exceeded = A. 61 x itfh^Mo or the dark matter mass in a cell is more than 2#idm> where mdm is the mass of the highest 
resolution dark matter particle. Cells are 610/i _1 pc across in comoving units at the maximum allowed refinement level. 

Five Milky Way size halos were simulated to z — 0. In order to isolate the effects of mass loss, simulations of each halo were 
run with mass loss turned on (ml) and mass loss turned off (nml). We required that the sample halos did not undergo major 
mergers after a lookback time of 5 Gyr and that sim ulations of the same halo with ml and nml did not differ significantly in 
accretion or timing of mergers(e.g.. lFrenk et alJ fT"999) as a result of the differing physics. Two halos met these criteria and they 
are presented here as MW1 and MW2. They have final total masses of M M wi = 1.6 X 10 12 r'M s and M MW 2 = 1.9 x 10 12 ft _1 M o 
at overdensity p = 200p cr , where p cl - is the critical density of the universe at the epoch of analysis. 

Subgrid Physics 

Star formation is modeled in our simulations with a simple density dependence and efficiency: 

. _ Pgas / Pgas \" 5 

P * ~ t. \0.01M o pc- 3 / 

, and we adopt t, = 4 Gyr, noting that the results and our conclusions are not very sensitive to the exact value of t». Stellar 
particles are allowed to form in regions with density above threshold value of n s f = 0.5 cm 3 and a temperature below the 
threshold value of T s f < 9000 K. Such a high temperature a nd low density thre sholds are reasonable choices for simulations of 
relatively low spatial resolution that we use in this study (see lSaitoh et al.l l2008. for a discussion of r», « s f, and T S f). 
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Each newly f orme d stellar particle is treated as an SSP with a stellar initial mass function (IMF) that is described by the 
iMiller & Scald d 19791) functional form, with stellar masses in the range 0.1 - 1OOM . All stars with M* > 8M immediately 
deposit both 2 x 10 51 ergs of thermal energy, accounting for the energy input by stellar winds and type II supernova (SN II), 
and a fraction f z = min(0.2, 2^Mt _ 0.06) of their mass as metals, into the surrounding gas, in order to approximate the results 
of iWooslev & Weaver! ([1995). The code also accounts for the SN la feedback, assuming a rate that slowly increases with time 
and broadly peaks at the population age of 1 Gyr. We assume that a fraction of 1.5 x 10~ 2 of mass in stars between 3 and 8 
M Q explodes as SN la over the entire population history and that each SN la dumps 2 x 10 51 erg of thermal energy and ejects 
1.3M of metals into the surrounding gas. For the assumed IMF, 75 SN II (instantly) and 11 SN la (over several billion years) 
are produced by a 1O 4 M stellar particle. This sort of model, without ad hoc kicks or delayed cooling, has been shown to have 
a negligible effect on the SFH and resulting galaxy properties compared to a model without feedback dKatzl [19921 [Tassis et al. 
2008; Sch aye et alJ l2010). Metallicity dependent equilibriu m cooling rates, as well as UV heating rates, due to the cosmological 
ionizing background, are tabulated using the Cloudy code (Ferland et al. 1998, v96b4). 

The stellar mass loss of each stellar particle is modeled using eq. Q] with Co = 0.05 and A = 5 Myr and about 40% of initial 
mass of stellar particle is lost over a Hubble time with these parameters. At each time step, gas mass lost by a stellar particle is 
added to gas of its parent cell along with its energy and momentum. The interpretation of our simulation results is not sensitive 
to the specific choice of Co and A. 

As a result of the known problems with simple galaxy formation simulations, our galaxies overcool and produce stellar fractions 
that are about 5 times too high compared to observational expectations. To account for this, we re-normalize all star formation 
histories of the simu lated galaxies w ith the constant ratio of the final stellar mass in the ml simulation to the Milky Way value of 
5.5 x 10 I0 M o re.g.. lFlvnn et"aT]l2006l) . 

Gas Reprocessing Results 

The star formation histories of our two simulated galaxies are plotted in Figure [8] There are no major mergers after z = 0.5, 
and these SFHs show few sustained bursts, making their growth largely consistent with expectations for the build up of the Milky 
Way (e.g.. lRocha-Pinto et al.l l2000). The SFHs show that z — SFRs in the runs with mass loss (ml plotted in red) are more than 
a factor of two higher than the SFRs in the corresponding simulations without mass loss (nml plotted in blue). This confirms the 
basic expectation that gas lost by stars can be a significant source of fuel for star formation at low redshift. 

We use the nml simulation to test whether simple assumptions about the fate of gas lost by stars can reproduce results of 
simulations with stellar mass loss. The GMLR is given by eq. [3] however, the gas returned to the ISM can accumulate and be 
consumed over a certain gas consumption time scale. The mass that accumulates in the ISM from stellar mass loss between time 
to and t is then 

m^(t ,t)= f {GMLR(f')-[SFR(f')-SFR nml (f')]}</r' (A1) 

J t 

To determine the SFR and hence a stellar mass loss rate, we assume a linear SFR scaling with gas density that has a gas con- 
sumption time scale r gc . Enhancement in star formation is therefore directly proportional to mass added to the star forming 
ISM, 

SFR(f) = SFR nml (f) + OTml( °' r) (A2) 

Stars formed from the gas lost by stars will themselves lose mass, an d this mass will, in turn form stars. In order to fully account 
for the additional star formation resulting from mass loss, equation IA2I needs to be iterated when starting from the nml case, 
replacing SFR nm i(f) with each subsequent SFR(f), until SFR(f) converges. 

The simulation's true star formation law is nonlinear in density, meaning that t sc effectively varies with density. The average 
gas densities of a star forming disk in simulations can vary from T gc = 500 Myr to a r gc of few billion years, with stars tending 
to form in lower density gas (longer r gc ) as accretion slows and the gas in their disks is depleted. However, this change in the 
efficiency of gas conversion occurs slowly and, furthermore, the effect of r gc is mitigated by the fact that a reduction in r gc results 
in a build up of m m \(0, t), which leads to more star formation. A constant r gc therefore mimics the relevant behavior. 

An alternative, simpler model, which we will call instant reprocessing 3 , can be constructed by taking t sc to 0. Then gas 
returned to ISM by stars is immediately consumed by star formation so only ongoing mass loss contributes to star formation and 
thus 

SFR(f) = SFR nml (f) + GMLR(r). (A3) 

We plot SFHs with both a r gc = 1 Gyr, and instantaneous reprocessing applied to the SFH from the nml simulation in Figure[8] 
The mapping between simulations is not exact. Gas returned during a merger may be returned to the diffuse halo and be temporar- 
ily lost to the star forming disk. Nonetheless, throughout most of its lifetime, and especially in late epochs when the galaxy is 
quiescently evolving, the simple assumption of instant reprocessing applied to the nml simulations generates an excellent match 
to the SFH in the ml simulations. 

Reprocessing models with longer gas consumption timescales result in gas accumulation and more gas reprocessing at late 
times, because they smooth the declining reprocessing rate curve. For the simulated SFHs, the effect of the r gc = 1 Gyr re- 
processing is a < 10% enhancement in SFRs at late epochs compared to the instantaneous reprocessing model. It is worth 



3 Instant reprocessing is distinct from the instant recycling approximation. Instant recycling approximates mass loss as instantaneous (e.g., Tinsley 1980), 
whereas instant reprocessing approximates the birth of stars from returned material as instantaneous. 
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Figure 8. Simulated SFHs that include all stars within 30/i 1 kpc of the galaxy and are normalized in such a way that the final stellar mass is equal to 5.5xl() 10 M o . 
The thin blue lines are SFHs for simulations that do not include mass loss from stars (nmt), while the thick red line is a fiducial simulation with mass loss included 
(ml). Black lines are gas recycling models, with gas consumption timescales r gc = OGyr (solid) and r gc = 1 Gyr (dotted), that attempt to map the nml SFH onto 
the ml SI'H. 

emphasizing that any feedback mechanism that slows the processing of gas to stars would have the effect of raising r gc and would 
therefore imply a higher rate of gas reprocessing at late epochs. 
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